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Anisotropic stress contributions to the gravitational field can arise from magnetic fields, colli¬ 
sionless relativistic particles, hydrodynamic shear viscosity, gravitational waves, skew axion fields 
in low-energy string cosmologies, or topological defects. We investigate the effects of such stresses 
on cosmological evolution, and in particular on the dissipation of shear anisotropy. We generalize 
some previous results that were given for homogeneous anisotropic universes, by including small 
inhomogeneity in the universe. This generalization is facilitated by a covariant approach. We find 
that anisotropic stress dominates the evolution of shear, slowing its decay. The effect is strongest 
in radiation-dominated universes, where there is slow logarithmic decay of shear. 


I. INTRODUCTION 

In this paper we consider the effects of anisotropic stresses on the late-time behavior of anisotropic and inhomo¬ 
geneous cosmological models which are close to isotropy and homogeneity. We use an approximation scheme that 
is similar to the velocity-dominated approximation but extend it to deal with the presence of anisotropic stresses. 
Past investigations of the late-time evolution of cosmological models in general relativity and in scalar-tensor grav¬ 
ity theories have been performed under the assumption that the matter source in the universe is composed only of 
isotropic perfect fluids. However, from the study of homogeneous anisotropic universes, it is known that the evolution 
of deviations from perfect isotropy is dominated by the distortions created by any anisotropic stresses. Examples of 
such stresses are provided by collisionless particles, gravitons, electric or magnetic fields, or by topological defects, as 
discussed in [jlj. 

In Q, a covariant and gauge-invariant approach was taken to the study of generic distortions to the Friedmann 
universe in the dust dominated era. The kinematic quantities which define deviations from perfect isotropy and 
homogeneity were related to the multipoles of microwave background temperature distribution. This enables limits to 
be placed on all forms of anisotropic and inhomogeneous distortions of the spacetime metric without making specific 
assumptions (like homogeneity) about the form of those distortions [||. 

In this paper we combine the approaches of jlj] and Q in order to provide a gauge-invariant and covariant description 
of the cosmological evolution of anisotropic and inhomogeneous universes containing a perfect fluid and anisotropic 
stress. This requires an extension of the covariant formalism developed earlier in Q to general anisotropic stresses and 
energy flux terms, and to the radiation era. The covariant approach is based on the full nonlinear equations, so that it 
naturally incorporates nonlinear effects. By the introduction of a self-consistent approximation scheme we determine 
the time-evolution of the anisotropic stresses and the ratio of the densities of the anisotropic to the isotropic fluid 
stresses. 

As in the case of spatially homogeneous universes [fj}|]j|] , the evolution of trace-free anisotropic stresses during the 
radiation era represents particular mathematical problems because of the nonlinear character of the anisotropic stress 
perturbations. If linearized about the exact Friedmann expansion this case would give rise to a zero eigenvalue. The 
dust evolution is determined in the linear approximation about the Friedmann solution but the decay of the shear 
anisotropy is considerably slowed by the presence of the anisotropic stresses. In each case there is an attractor solution 
relating the distortion in the expansion anisotropy to the fractional density in anisotropic stress. Here we find that 
similar behavior persists in the presence of small inhonrogeneity. 


II. COVARIANT LAGRANGIAN DYNAMICS 

The source of the gravitational field is a non-interacting mixture of perfect and anisotropic fluids, the former 
typically representing isotropic radiation or cold matter (dust), while the latter encompasses an ultra-relativistic 
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fluid, collisionless or viscous, a magnetic field, and other possibilities. The total energy density is 

Ptot — Pperf d“ Paniso = P T P ; 


and the total isotropic pressure is 


Ptot = (7 ^ !)P + \P > 


where 7 is constant, with 7 = | for radiation and 7=1 for dust. The perfect fluid is dynamically isotropic about a 
unique 4-velocity u a , where u a u a = — 1. Relative to u a , the energy flux of the anisotropic fluid is q a and its anisotropic 
stress is 7 r a b, so that the total energy-momentum tensor is 

Tab = {p + p)u a Ub + [(7 - T)p + \p]h ab + q a u b + u a q b + 7r ah , (1) 

where h ab = 9ab + u a u b is the projection tensor into the instantaneous rest spaces of comoving observers, with g ab the 
spacetime metric. 

The anisotropic dynamic quantities q a and T: ab can in general encode either relative velocity effects or dissipative 
effects, or both. The total energy flux in general comprises a heat flux Q a and a particle diffusion flux (momentum 
density) [§ 

Qa — Qa ~ 1“ 3 f^Va 5 

where v a is the velocity relative to u a . In the case of a magnetic field, ir ab arises from the anisotropic contribution 
of the field to pressure. For a fluid in the hydrodynamic regime, ir ab is driven by shear viscosity. For a collisionless 
gas, ir ab arises from the covariant quadrupole moment of the distribution function. The full anisotropy information is 
encoded not only in q a and n ab , but also in the octopole and higher covariant multipoles of the distribution function 

i- _ 

The covariant Lagrangian approach |7J] is based on a physical 4-velocity u a , with all physical quantities represented 
by tensors in the comoving instantaneous rest space. This involves a 1 + 3 splitting via u a and h ab . The covariant 
time and spatial derivatives areQ S a "' b -.- = u c V c S a '" b - - and T) c S a ' b ... = hjh a d ■ ■ ■ h b e ■ ■ ■ VfS d '" e .... (Note that 
D c h ab = 0.) As part of a streamlining and development of the covariant formalism, we define a covariant spatial 
divergence and curl of vectors and rank -2 tensors |)J : 

div V = D q R q , (div S) a = D b S ab , ( 2 ) 

curl Va = £abcT) b V c , curl Sab = £cd( a T> c S b ) d , (3) 

where £ a bc = Pabcd,u d is the covariant spatial alternating tensor, projected from the spacetime alternating tensor 
Pabcd = — (It satisfies £ abc £ def = 3 lh[ a d h b e h c ]f.) These operators are irreducible parts of the 

spatial derivative, making transparent the electromagnetic structure of the general relativistic equations, and they obey 
covariant identities which greatly facilitate manipulation of the equations to investigate perturbations, integrability, 
wave equations, and so on (see | 0 ). 

A projected vector satisfies V a = V( a ) = h a bV b . The covariant irreducible decomposition of a spatial rank-2 tensor 
is S ab = \S cd h cd h ab + S( ab ) +£ abc S c , where the PSTF part is S( ab ) = h( a c h b ) d S c d - \S cd h cd h abl and the skew part is 
the dual of the spatial vector S a = \£ a be s [6c] . 

The kinematic quantities arise from decomposing V b u a : 

^bUa = DfrUa A a 1i b , (4) 

1T &ab T £abc^ ; (5) 

where A a = u a = A^ is the acceleration, 0 = D a w Q is the volume expansion rate, a ab = D( a Ufr) is the shear, and 
oj a = —icurl u a is the vorticity. The Weyl tensor, representing the locally free gravitational field 0, is decomposed 
into gravito-electric and magnetic fields: 


Tab — C acbd U U — ^(ab) ? T[ ab — 2 £acdC be'U — -^(afo) • 


( 6 ) 


1 In Q, D a is denoted V a , and in [^] it is < - 3 * V„. 
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Then the background Friedmann model is covariantly characterized as the limiting case 


dynamics: D a p = 0 = D Q p, , q a = 0 , Tr ab = 0 , (7) 

kinematics: D Q 0 = 0 , A a = 0 = u> a , a ab = 0 , ( 8 ) 

free field: E ab = 0 = H ab . (9) 

In covariant and gauge-invariant perturbation theory j|], the above quantities are non-zero but small, so that the 
universe is almost isotropic and homogeneous. In spatially homogeneous anisotropic models (Bianchi and Kantowski- 
Sachs models), out of the quantities in equations ( 0 )-®, q a , t Tab, cr ab , E ab and H ab may be nonzero. 

The field equations R ab = T ab — ~T c c g ab , with T ab given by Eq. ( 0 ), are incorporated into the Ricci identity 
V[ a Vfc]W c = \R a bcd,u d and the Bianchi identities V d C a w = V[ Q (— R b j c -f \Rg b ] c ) 1 which are then covariantly decom¬ 
posed into propagation and constraint equations as in [jjj. The equations are significantly simplified and made clearer 
in the developed version of the formalism |J. This gives the following (exact, nonlinear) equations: 


Propagation: 


p + 1 Qp = 0, (10) 

/i + |0/i + D a < 7 a = —cr ab TT ab — 2A a q a , (11) 

< 7 <a> + I©< 7 a + §Ma + + (div 7 t) q = -(<T ab + £ abc UJ C )q b - A b TT ab , (12) 

0 + ^0 2 + \ [(3y — 2 )p + 2p] — D a A a = —a ab cr ab + 2 oj a u) a + A a A a , (13) 

W( a ) + §0w„ + |curlA a = a ab io b , (14) 

d (ab) + 10C a b + E ab — \TT ab — D( a A b ) = —Cr c ( a (Tb) C — W( a W b ) + A( a A b ) , (15) 

E(ab) + QE ab - curl H ab + \n( ab ) + ^07 Tab 

+ lilP + | P) a ab + {a q b) = Vc{a [3E b) c - ±7r h) c ] - u c £ cd ( a [E b) d + ^TT b ) d ] 

— A( a qb) + 2 A c £ cd ( a H b ' > d , (16) 

R(ab) T 077ab T CUrl E ab 5 curl7r ab — 3(T c( a E b ) ^ ^cd(aRb) 2A £cd(aE b ) 

2^( a Qb) 4 “ 2 ^ (a^b)cd,q • ( 17 ) 


Constraint: 

"1pA a + (7 - l)D a p = 0 , (18) 

§D a 0 - (div cr) a + curl u) a - q a = 2 £ abc u b A c , (19) 

&rvu] = A a uj a , (20) 

Rab CUll (Job Y)^ a CJ b ^ — 2 A^ a iO b ) j (21) 

(div E) a - ±D a (p + n) + lD b 7r a6 + ±0g a = £ a bcP b dH cd ~ 3H ab LO b + \cr ab q b + § £ a b c u b q c , (22) 

(div H) a - (jp + | p)u a + \ curl q a = - £ a bcV b d [E cd + §7r cd ] + [3 E ab - ±7r ah ] uj b . (23) 


We emphasise that these are the exact equations of the full nonlinear theory, before any approximations are made. 
They generalize the equations given for the case of irrotational dust in |j| and perfect fluid in [|o| , and the linearized 
dissipative fluid equations given in Q. Linearization about a limiting background Friedmann model reduces all of 
the right hand sides to zero, in view of equations ©-(!)• (Linearization also allows us to remove the angled brackets 
from the indices of time derivatives.) 

We note that if A a = 0, then Eq. (|l^) implies D a p = 0 unless 7 = 1. It then follows from the energy conservation 
equation © (recall that the two fluids are non-interacting) and the identity 

D a <fi = (D a (j>Y + 10D a <l> + cr a b V b 4> + £ab c u b D c </> - A a (j> - u a A b D b cj >, 

that D a 0 = 0 also. Furthermore, the identity 

curlD a <£ = - 2 (jxjj a , 
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shows that Lo a = 0. In summary: if the perfect fluid is not dust, then vanishing acceleration (A a = 0) implies vanishing 
vorticity (uj a = 0) and spatial gradients of energy density and expansion rate (D a p = 0 = D a @). 

The anisotropic stress ^x ab may be subject to an evolution equation, i.e. it may be self-consistently determined via 
kinetic theory || or thermodynamics |j,[ll| or some other dynamics |l],|l2|,|3|. For example, in the case of a collisionless 
gas, the covariant quadrupole of the Boltzmann equation is a propagation equation for Tt ab & 

K(ab) + f ©7 Tab + JgPCTab + 2D {a q b) + (div£) ab = 

- 2 A {a q b) - |CT c ( a 7r b)c + 2u d £ cd{a Tx b) c + a cd Xabcd , (24) 

where £ abc and Xabcd are the octopole and hexadecapole. Another example is a dissipative fluid, for which the causal 
thermodynamics of Israel and Stewart gives Q 

TTT(afc) + 7 Tab = ~2T]a a b , (25) 

where 77 is the shear viscosity and r is the relaxation time-scale, with r = 0 in the non-causal Navier-Stokes-Fourier 
theory. The heat flux in such a fluid obeys the causal transport equation^ 

T Q(a) + Qa = — k (D 0 T + TA a ) , (26) 

where k is the thermal conductivity and T is the temperature. 

In this paper, we will not impose evolution equations such as (|24|) (0), since our main interest is in the general 
dynamics of anisotropic stress, rather than particular physical models of stress. For this purpose, we will introduce 
in the next section a general expression for 7r ab . 


III. APPROXIMATIONS AND SOLUTIONS 


The equations of the previous section provide in principle a complete, fully covariant and nonlinear description of 
the effects of anisotropic stresses on cosmological evolution, purely in terms of quantities with a direct and transparent 
physical interpretation. However, the complexity of the exact nonlinear system effectively rules out such a description, 
and we must resort to approximations. The fundamental approximation is that the universe is close to a Friedmann 
model. This in fact follows after last scattering as a theorem |14| , using the observed small anisotropy in the microwave 
background temperature, and assuming a weak Copernican principle (which is falsifiable via the Sunyaev-Zeldovich 
effect |j)). Before last scattering, in the radiation era, it is reasonable to make the assumption because rapid collisions 
should tend to damp out perturbations in the dominant component, leaving anisotropic stresses of the decoupled 
component as a perturbation. 

As pointed out above, it is not sufficient to linearize the equations about the Friedmann limit, since during radiation 
domination the anisotropic stress perturbations have an essential nonlinear character, which is known in homogeneous 
anisotropic models and will be confirmed to extend to the small-inhomogeneity case below. A straightforward 
linearization would reduce all the right hand sides of the propagation and constraint equations (10)-(|23|) to zero. In 


order to incorporate the lowest order effects of anisotropic stress and energy flux in the radiation-dominated era, we 
must retain on the right hand sides those terms which couple Tr ab or q a to the shear cr ab . 

Even when we neglect all second-order terms apart from those coupling shear to anisotropic stress and energy 
flux, the equations remain formidable. Since our main interest is in how the late-time behavior of shear is affected 
by anisotropic stresses on large scales, we impose the so-called velocity-dominated approximation, i.e. that spatial 
derivatives are neglected with respect to time derivatives. This means that in the propagation equations 
we neglect spatial derivative terms relative to the other terms. 

We normalize the key quantities so as to obtain a set of covariant dimensionless quantities. These are chosen as 
follows: 


* D aP * Qa * 

Pa 11 5 Qa ’ n ab 

Hp p 

A* — — ui* - — 
A a - H^a- H 


* M 
P = ~ 

P 

H 2 

p* _ E ab _ H ab 

^ab — ’ n ab — ’ 


7Tab 


a ab ~ 


®ab 


2 This neglects any coupling of heat flux and viscosity, and assumes the relaxation time-scales for heat conduction and shear 
viscosity are equal. For the general case, see j(|. 
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where H = a/a is the background Hubble rate. We also need the background field equations 


p = 3H 2 (1 + K ), H = ~\H 2 [ 3 7 + (3 7 - 2 )K \, 

where K = JC/(aH) 2 and K. = 0, ±1 is the curvature index. The above dimensionless quantities are small, and vanish 
in the background. Thus they are gauge-invariant as well as covariant Jsj. 

Using these quantities, applying the velocity-dominated approximation, and retaining terms that couple shear to 


anisotropic stresses, the propagation equations imply 

A* = -(4-3 7 )Hp*~Ha* ab n* ab , (27) 

q: = -(4~^)Hq:-Ha: b q* b , (28) 

^ = -H 4 - 3 7 + (2-3 7 )i^:, (29) 

Kb = -H 4 - 3 7 + (2 - 3 j)K\ Ha* b - HE* ab + f (1 + K)Hn* ab , (30) 

K b = - [3(1 - 7) + (2 - 3 7 )K] HE* b - | 7 (1 + K)Ha* ab - |(1 + K)n* ab 

+ |(3 7 - 1)(1 + K)Hn: b - |(1 + K)Hn* c{a K } K (31) 

H* ab = - [3(1 - 7 ) + (2 - 37 )K\ HH* ab + |(1 + K)Ha* c(a e b) cd q* d . (32) 


Eq. ( p 8 | ) is a propagation equation for the energy flux. There is no propagation equation for the anisotropic stress, 
and we need a general prescription which allows us to investigate qualitative features of the role of anisotropic stress 
on shear dissipation. Following previous studies in homogeneous universes jl||],|l5|], we assume that the anisotropic 
stresses are of the form 


'ab 


= p*\ 


ab 5 


^ab — 0 ? 


(33) 


where A a b is a tensor whose components are slowly-varying functions of the space coordinates x relative to the time 
variation of p*. 

We shall further assume that the universe is close to a flat Friedmann model, so that K = 0. The self-consistent 
approximation for such an analysis is that the mean Hubble rate and the perfect fluid density evolve, to leading order, 
as in the Friedmann universe, so that 


H = 


3 jt 

Then the non-linearly coupled evolution equations 

( 4 - 3 7 )' 


P = 


3 7 2 f 2 
for p* and 


for a* b become 


2(4 - 3 7 ) 


a ab 


3 jt 


a ab 


3 jt 

(— 

V 3 7 * 


p = 


Kb = 


_3 jt 

K_ 

It 


ab„* 
ab 


-A a. 


a b • 


(34) 


(35) 

(36) 


First, we note that a qualitatively different situation arises according as 7 =|or 7 ^|. If | < 7 < |, then 
the right-hand side of Eq. (^) can be neglected, and the evolution of p* is determined at linear order in p* by the 
solution 


p* (t, x) = M(x)t n 


where 


n = 


2(4 — 3 7 ) 
3y 


(37) 


where M is a slowly varying function of x. This expression also holds formally when 7 > |, but, since n > 0 in that 
case, the effects of the anisotropic stresses grow in time and the assumption ( |34| ) that the mean expansion is close to 
the Friedmann model, must eventually break down. Our analysis does not apply to this case. When 0 < 7 < |, the 
anisotropic stresses do not dominate over the simple cr oc t~ 2 ^ mode that arises when p* = 0. When | < 7 < |, 
p* decreases with time as t —► 00, and the approximation becomes increasingly better. However, the fall in cr* is 
always at a slower rate than occurs in the absence of anisotropic stresses (p* = 0). Note that when 7 = 1, we have 
p* oc f -2 / 3 . 

However, when 7 = |, the situation changes. The evolution of p* is then determined by the nonlinear terms 
on the right-hand side of (|35|). In this situation the universe evolves towards an attractor in which a* b varies only 
logarithmically. 
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For general 7 < Eq. (p 6 |) has relaxation form, and we seek asympotic solutions with a* b —> 0. This gives 


Hence for 7 < | 


£4 = §^A Q6 + i(3 7 -4K 6 . 


Kb °c <6 oc /X* OC r 


( 38 ) 


(39) 


where the factors of proportionality are slowly-varying functions of the space coordinates. 

4 
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Specializing to the radiative case 7 = ^ again, the evolution of the gravito-electric field is determined by Eq. 


K b = (^W- ()W+ (| 


f - (| 


) . ,* \ c _* 

A 4 A (a a b)c 


(40) 


Combining this with Eq. (^ 8 |), neglecting all but leading order terms in the shear, we have in this approximation 

Kb = 3KX ab - (41) 

Then, together with Eq. (p5|), this implies 


_ _ ° ,.*2 \ \ab 

M ■ 


and therefore, 


^ ^ 3X a b{x)X ab (x) [IV(x) + lnt] ' ^ ^ 

where N(x) is an arbitrary slowly-varying function of integration. Thus, a* b is determined in this approximation from 
Eq. (f4i|) and E* b is determined from Eq. (|3|). Asymptotically, as t —> oo, we have 


Kb oc Kb °c Kb °c f cx jE ) ( 43 ) 

c Tab OC -E , E ab OC 7 Tab OC d ^ . (44) 

t In t V In t 

Moreover, in this approximation, using equations (|3^) and ©> we have that at large t, 

E a b = \^ab- (45) 

Note that the shear evolution found here is dictated by the anisotropic stresses. If we were to take TT a b = 0, as is 
usually done in analyses using the velocity-dominated approximation, then A a b = 0 and we do not have cr* b —> (lnt ) -1 
asymptotically. Instead, we recover the familiar Bianchi type I mode with a a b oc a -3 cx t -3 / 2 , where a is the 
geometric-mean expansion scale factor (H = a/a). 

The propagation equation (|2|) for energy flux becomes, to leading order, 

fa = - Ha* ab q * b . 

Hence, in the 7 = | case we have, as t —> oo, 


It follows that 


1 

( X cd X cd ) tint 


A a b q* b ■ 


* 

Qa OC 


1 

lnt 


O q a OC 


1 

t 2 In t ’ 


( 46 ) 


so that the energy flux has the same asymptotic behaviour as the anisotropic stress. 
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IV. CONCLUSIONS 


In this paper we have set up an efficient covariant formalism for studying the asymptotic behaviour of cosmological 
models. This formalism allows us to delineate the contributions from linear and non-linear terms in the evolution. 
We have investigated the influence of trace-free anisotropic stresses on the late-time evolution of universes containing 
small anisotropies and inhomogeneities. These stresses can be contributed by collisionless particles, magnetic fields, 
or topological defects and are found to dominate the evolution of the anisotropy. In the radiation-dominated case 
the contribution of trace-free anisotropic stress tensors is found to be mathematically subtle. The perturbations to 
isotropy appear at second-order in the shear and lead to a slow logarithmic decay of the shear anisotropy. In the dust 
era the fall-off is as a power-law in time but is again slower than it would be in the absence of anisotropic stresses. 
These results generalize those found earlier for homogeneous anisotropic models by Barrow jj]|| (see also [|l5j]). They 
extend the analysis of Maartens et al || to the case where the background model model contains fluid with non-zero 
pressure. These analyses show us that asymptotic studies of the late-time evolution of expanding universes fl(if| need 
to include the effects of anisotropic stresses in order to discover the leading-order form of the metric. This will also 
apply to studies of the late-time evolution of cosmological models in other theories of gravity, for example those in 
scalar-tensor gravity theories performed by Cromer et al |l7j. 
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